Entanglement generation by collisions of quantum solitons in the Born's 

approximation 
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We present analytic expressions describing generation of the entanglement in collisions of initially 
uncorrelated quantum solitons. The results, obtained by means of the Born's approximation (for fast 
solitons), are valid for both integrable and non-integrable quasi-one-dimensional systems supporting 
soliton states. 
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I. INTRODUCTION 

Entanglement is a fundamental property of bipartite 
quantum systems Apart form being a major resource 
for quantum-information techniques, entanglement ex- 
hibits itself perhaps in the most spectacular form in the 
breakdown of Bell's inequalities and Einstein-Podolsky- 
Rosen "paradox" . In experiments, strongly entangled 
states, and in particular those exhibiting "nonlocality" , 
are typically created with microscopic particles produced 
by the same source, or interacting prior to the detection, 
such as pairs of photons 0, ^ or ions Q . In their orig- 
inal paper Einstein, Podolsky and Rosen considered 
two particles whose ordinary degrees of freedom, center- 
of-mass positions and relative momenta, were correlated. 
This situation is realized with the entanglement of Gaus- 
sian states of light [6] or atomic ensembles Q . 

Particularly interesting is the possibility to generate 
entanglement between macroscopic (or mesoscopic) ob- 
jects that may be transmitted over long distances. In 
this connection, solitons, i.e., stable solitary waves that 
propagate without distortion may be considered as 
possible robust quantum-information carriers. In partic- 
ular, the solitons as collective excitations in nonlinear 
media (unlike ions or other material objects) may be cre- 
ated with a desirable shape and effective mass, and ad- 
mit a much greater degree of control by means of various 
"management" techniques 

The objective of the present work is to analyze the 
possibility of the creation of entanglement between soli- 
tons by means of collisions between them. To the best 
of our knowledge, this problem was not considered be- 
fore (except for a preliminary version of the present work 
p^). However, a very recent work ^13j has presented 
a rigorous proof of the generation of entanglement be- 
tween constituent solitons in oscillating two- and three- 
soliton bound states, in the integrable model based on 
the NLS (nonlinear Schrodinger) equation with the cu- 
bic nonlinearity. It is relevant to note that multi-soliton 
solutions of the integrable NLS equations are actually un- 
stable against slow separation, hence the entanglement, 
generated by the interactions between the solitons while 



they remained bound, may be kept after the separation. 

Natural candidates for the study of the collision- 
induced entanglement are matter- wave solitons, that can 
be built of ultracold bosonic atoms, which form the Bose- 
Einstein condensate (BEC) Q, and photonic solitons 
in nonlinear optical waveguides d, [lOl . The studies of 
BECs has led to the creation of dark [l4l [Tst and bright 
solitons in trapped Bose-Einstein condensates 
(with repulsive and attractive interactions, respectively). 
Recently, the observation of stable bright-dark soliton 
pairs has been reported too (21i |. In addition, similar- 
ities to nonlinear optics Q have triggered the interest 
in theoretical studies of discrete (lattice) BEC solitons 
[12, , and have also led to the seminal observation of 
gap solitons, i.e., robust localized matter- wave packets 
supported by the interplay of repulsive interactions and 
an effective negative mass of collective excitations in the 
condensate, induced by the periodic optical-lattice poten- 
tial [205. The analysis reported below suggests that the 
matter- wave solitons in BEC with attractive inter-atomic 
interactions have the best potential for the generation of 
the entanglement through collisions; in particular, the ve- 
locity of the moving solitons, which is the crucial parame- 
ter, which determines the collision- induced entanglement 
(see below), can be easily controlled for BEC solitons. 

While most of the previous studies of matter- wave soli- 
tons were concentrated on their classical and mean-field 
aspects, more recently considerable interest has been de- 
voted to the role of thermal noise 23, 25] and quan- 
tum fluctuations. The latter may cause filling up of the 
dark-soliton's core through the quantum depletion pro- 
cess, as was predicted using the Bogoliubov-de Gennes 
(BdG) equations [1^. Making use of the discrete NLS 
equation and the time-evolving block-decimation algo- 
rithm '275 in the framework of the Bose-Hubbard model, 
it was confirmed that quantum effects lead to the filling 
in of dark BEC solitons, and it has been demonstrated 
that collisions between them become inelastic [28j. The 
BdG approach and its generalizations were also employed 
to study the solitons' stability [29, 30, 31] and excitations 
caused by opening of the trap [3^. In the latter work, 
exact Lieb-Liniger solutions were also used, to calculate 
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internal correlation functions of positions of the particles. 
A noisy version of standing bright solitons was studied by 
means of the exact diagonalization and quantum Monte- 
Carlo method [1^. 

Being typically associated with integrable ID classical 
models, stable solitons are known too in non-integrable 
systems, including multidimensional ones. One of re- 
markable properties of solitons is their stability to pertur- 
bations. The stability can be extended to quantum set- 
tings, in which the mean- field description admits quan- 
tum solitons in the semi-classical form. As mentioned 
above, the robustness of mean-field solitons against quan- 
tum fluctuations (including finite-temperature effects) 
was studied, using the time-dependent Hartree-Fock- 
Bogoliubov equations, for both ordinary matter- wave 
solitons, supported by attractive interactions between 
atoms [13], and for gap solitons [Slj. Although under 
extreme conditions the quantum fluctuations may split 
a mean- field soliton (s^l, it has been concluded that, 
in a broad range of parameters relevant to the experi- 
ments, the matter-wave solitary waves predicted by the 
mean- field description (i.e., found as stable localized so- 
lutions to the Gross-Pitaevskii equation, GPE are 
completely robust objects - in fact, in perfect agreement 
with the experimental observations of these solitons p7| - 
[20| . As concerns the relation between the mean- field 
and quantum descriptions of localized objects in BEG, 
it is relevant to mention that an alternative derivation 
of the GPE from a consistent many-body quantum the- 
ory, based on the variational approach in the multi- 
configurational space, was recently presented in Refs. 

In this work, we analyze the collision-induced gener- 
ation of the entanglement in pairs of fast solitons for a 
class of equations of the NLS type by means of the Born 
approximation, which is valid when the kinetic energy 
of the moving objects (solitons) is much larger than the 
potential of the interaction between them, hence the in- 
teraction may be treated as a perturbation [3^ (there- 
fore, this is the case opposite to that studied in Ref. p^ . 
where the entanglement was considered between bound 
solitons created with zero relative velocity). We aim to 
present simple analytical expressions for the collision- 
generated entanglement of two quantum solitons, which 
are valid for generic quasi-lD systems - integrable or not 
- that admit soliton solutions. Another approach is possi- 
ble for nearly quiescent solitons, when the entanglement- 
generating perturbation is the weak interaction between 
them, assuming that the distance between the solitons is 
large enough. Recently, this approach was developed in 
Ref. [S^l for kinks in the sine-Gordon equation. In that 
model, the kinks were maintained in the quiescent state 
by boundary conditions, as the model was defined in a 
domain of a finite size. 

The rest of the paper is organized as follows. In Sec- 
tion II, we introduce the model, starting with the known 
system of classical equations of motion for a pair of well- 
separated solitons [10], and then proceed to the respec- 



tive quantum system. We confine the analysis to the ba- 
sic case of the pair of symmetric solitons with equal am- 
plitudes. Mismatch between the amp litudes makes the 
interaction effectively incoherent [3a], thus suppressing 
the entanglement generation (we estimate the size of the 
mismatch up to which it may be neglected). The initial 
quantum state, corresponding to far separated solitons, is 
taken as a product of two independent wave packets. In 
Section III, the calculation of the correction to this fac- 
torized state, generated by the collision between the fast 
solitons, is performed in an analytical form, by means of 
the above-mentioned Born approximation. The collision- 
induced correction to the wave function features explicit 
entanglement in terms of two relevant degrees of freedom, 
viz., the distance between the solitons and their relative 
phase, r and x- The results are reported for the initial 
condition of two types: a more sophisticated one, with 
the Gaussian localization of x around a definite value, 
and also for a simple phase-independent initial distribu- 
tion. Both types of the initial quantum states may be 
realized in the experiment, under different specific con- 
ditions. The paper is concluded by Section IV, where, 
in particular, we discuss the robustness of the predicted 
entanglement against external noise, and possible exten- 
sions of the work. 



II. THE MODEL 

A. The classical soliton pair 

For the condensate with attractive interactions be- 
tween atoms, the scaled form of the GPE in the free ID 
space, which describes the BEG in the mean-field approx- 
imation, is tantamount to the integrable NLS equation 
[H, i.e., 

iut + {l/2)u^^ + \u\'^u = 0. (1) 

The commonly known soliton solution to Eq. ([T|) is 

Usoi = ?/sech [i]{{x — S,(t))] exp [i(l)(t) + icx] , (2) 

where rj is the soliton's amplitude, d^/dt = c its velocity, 
and 

- d(j>/dt = (c^ - jf) /2 (3) 

the intrinsic frequency. 

More realistic forms of the GPE in ID include vari- 
ous terms which break the integrability of the NLS equa- 
tion. In particular, the full equation must include the 
axial trapping potential, but, as concerns interactions 
between solitons, the axial potential is not a crucially 
important factor, according to the available experimen- 
tal results [13] and theoretical analysis [H, |4l[ . Another 
physically relevant feature that may affect soliton-soliton 
collisions is a quintic nonlinear term. It may account for 
three-body collisions in the condensate, provided that 
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they are lossless [42|, but a more general (in fact, uni- 
versal) source of the quintic term is the deviation of the 
condensate loaded into a cigar-shaped trap from the one- 
dimensionality. This term can be derived by means of a 
perturbative analysis ^43] , or by the expansion of a more 
general equation, that takes into regard the underlying 
three-dimensionality via the nonpolynomial nonlinearity 
p^ l- The universal quintic term always corresponds to 
self- attraction (irrespective of the sign of the binary in- 
teractions between atoms), with the coefficient in front of 
it proportional to the transverse-confinement frequency 
and square of the coUisional scattering length. By means 
of straightforward rescalings, the equation with the com- 
bination of attractive cubic and quintic terms can be cast 
in the form that contains no free parameters, 



iut + {l/2)uxx + \uf u + \u\'^u 



0. 



(4) 



Although the cubic-quintic (CQ) NLS equation (|4]) is not 
integrable, its exact soliton solutions are well known [ist . 
the entire family being stable [46| : 



y/2r] exp [i0(i) -I- icx] 



^^1 + (8/3) ry2 cosh pry {x - 



(5) 



where 0(i) and £,{t) have the same meaning as in Eq. ([2]). 
The CQ model based on equatio n ffl finds other physical 
realizations in nonlinear optics [47[ . 

The analysis which treats the interaction between well 
separated quasi-classical (mean-field) solitons ([2|) with 
common amplitude 77 as a perturbation leads to the the 
conclusion that they may be considered as a pair of quasi- 
particles with an effective Hamiltonian [lo| . 



Hu 



(1/2) (p' 



a2 



X 



cosx, 



(6) 



where the overdot stands for the differentiation with re- 
spect to rescaled time r = 2y/2r]^t, the normalized dis- 
tance between the solitons and center-of-mass coordinate 
are defined as 



= ^(a-6), P = ??(Cl+6): 



(7) 



while the relative and overall phases of the soliton pair 
are 



X = 01 



+ 02 - V t. 



(8) 



A noteworthy fact, which is obvious in expression ([5]), is 
that the effective mass corresponding to the phase de- 
grees of freedom in the two-soliton set is negative [48] . 

Actually, Hamiltonian ^ is universal, applying to any 
model which supports exponentially localized solitons 
with an intrinsic phase degree of freedom px| . In partic- 
ular, this Hamiltonian governs the interaction between 
solitons ^ of the NLS equation with the CQ nonlinear- 
ity, Eq. ([3]) , up to a proper rescaling of coefficients ^] . 



B. Quantization 

Proceeding to the quantum version of the model con- 
sidered above, we treat, as usual [ioj, each soliton as a 
quantum particle with two degrees of freedom, the posi- 
tion and phase (in the experimentally relevant situation, 
effects of quantum fiuctuations around the quasi-classical 
shape of the solitons may be negligible for matter-wave 
solitons, as discussed above). Thus, the quantum coun- 
terpart of classical Hamiltonian [5] gives rise to the follow- 
ing linear Schrodinger equation for the wave function of 
the soliton pair, vI/('°*), which depends on the total set of 
four degrees of freedom describing the pair: 



(tot) 



'(1/2) r^'(*°*) - -)- ^f^*"') - ^1*°*'' 

V^/^/ ^ rr XX PP 06 

-e e-l''! (cosx)*''°*\ 



(9) 



where T = Ht and e = l/ti^, with h the renormal- 
ized Planck's constant [measured in scaled units in which 
Hamiltonian ^ was written] . The reduced form of Eq. 
© for the wave function which depends only on the rel- 
ative variables, r and x, is 



i^T = - 



^(*r^-*xx)+£e I'^lcosx 



(10) 



Our objective is to analyze collisions between rapidly 
moving quantum solitons, with large relative momentum 
Kf). The initial state is taken as a naturally expected 
non-cntangied factorized one, centered around definite 
initial values of the dynamical variables, ±^0 and ^Xo- 

(6,6,01,02) = 



exp 



X exp 



Va (6 - ^of (01 - Xo)^ 



2S2 

„2 . - n2 



2$2 



vl (6 + Co)- (02 + xo)^ 



2S' 



"2 



2$2 



(11) 



This state assumes equal widths of the wave packets for 
the two solitons, S and <i>, with the pair's center of mass 
set at X = 0, and initial separation 2^ = £,i{T = 0) — 
6(7" — 0). The mean value of the initial overall phase 
is also fixed to be zero, while the initial phase difference 
between the solitons is 2xo = 0i(T = 0) — 02 (r = 0). 
Finally, initial state pT|) can be written in terms of the 
variables defined in Eqs. ([7]) and ([5]), 



X exp 



(6,6,01,02) = 



4S2 



4$2 



(12) 



where and = 2r/6- 

To conclude the formulation of the model, it is rele- 
vant to consider the possibility of the excitation of an 
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intrinsic mode (IM) in the colliding solitons (if the IM 
exists). A known principle is that solitons of integrable 
equations do not support IMs, but a nonintegrable model 
may feature an IM. In particular, exactly one IM is exists 
in excited states of solitons ^ of the CQ NLS equation, 
with the self-focusing sign of both nonlinear terms [i^ . 
A possibility of using solitons' IMs as carriers of quan- 
tum information was proposed in Ref. [37]. However, 
in the case of the collision between fast solitons, which 
is considered below, the excitation of the IM, as well as 
generation of nonsoliton modes ("radiation"), may be ne- 
glected, in the lowest approximation, simply because the 
intensity of these effects is inversely proportional to the 
square of the coUision velocity [lol |. 



III. GENERATION OF ENTANGLEMENT BY 
COLLISIONS BETWEEN SOLITONS 

A. The Born's approximation 

Scattering solutions to reduced equation (|10[) are 
generated by the incident wave, ^'(r —^ oo) = 
exp ^iKr + inx — i [K^ — k^) T/2] , where k is an inte- 
ger. The full solution is sought for as 

* (T, r,x)^V (r, x) exp [-i (K^ - k^) T/2] , (13) 

with the stationary part of the wave functions obeying 
the following equation, 



Vrr - V, 



XX 



2£ e-l''lcosx+ V^O. (14) 



As said above, our basic assumption is that we consider 
the collision between fast solitons, i.e., is a large pa- 
rameter in comparison with e {k^ may be large too). In 
other words, the kinetic energy of the relative motion is 
much larger than the potential of the soliton-soliton in- 
teraction. In fact, the two-particle Hamiltonian ^ can 
be used for the description of collisions between solitons 
only in this case; otherwise, one cannot neglect deforma- 
tion of the solitons in the course of the collision, as well 
as the excitation of the IM, if it exists in the solitons, 
and the generation of the radiation modes. 

We apply the Born's approximation (sH], looking for a 
solution to Eq. as 

V{r, x) = [1 + <yJ (r, x)] exp {iKr + inx) , (15) 

where perturbation which is assumed to be a slowly 
varying function of r in comparison with exp (iKr), obeys 
the simplified equation: —iKipr + {l/2)(p^^ + inip^ = 
£ e~l''lcosx- Solutions to this equation satisfying the 
necessary boundary condition, ip{r = — oo) = 0, can be 
readily found in an analytical form. The outcome of the 
collision is determined by the asymptotic form of the so- 
lution, which is found to be 



f{r +00, x) — (is/K) 

X g»Xgi(K-l/2)r/X ^ g-iXg-«(K+l/2)r/x' 



(16) 



B. Analysis of the colHsion-induced entanglement 

In the case of the fast collision, the overall variables, 
p and 9, may be treated as "frozen" ones in the wave 
packet generated by initial state ([T^ . Then, it is natural 
to decompose the initial state over the set of plane waves 
with respect to variables r and x' 



4-^'^ (^1,6,01,02) = exp 



452 4$2 ) 



(S$/7r) / dK d«;e'-^('-''«)+''(x-xo) 

J — 00 J — 00 

xexp\-E^ (K - Kof ^^K^] . (17) 



Note that the Gaussian distribution of angular wavenum- 
ber K, which must be integer, is valid for |k| 3> 1, i.e., 
$ <C 27r [recall that $ is the width of the initial distribu- 
tion of the phase variables introduced in Eqs. pT|) and 

mi 

Next, recombining wave packet pT|) with the collision- 
induced perturbation of the wave function, as per Eqs. 
(fT5|l and (fT6|) . and again making use of the fact that Kq 
is large, we arrive at an expression for the net change 
of the wave function which is generated by the fast col- 
lision between the two solitons, in the first order of the 
perturbation theory: 



<5* (ei,6,'/'i, 



*2)«exp(-^- — 



Ko (4$4 + T2) 
e^*^ exp 



1/4 



exp 



<i>^ {X~X0± r/KoY 

4$4 ^ y2 



(18) 



with phases shifts 



tan" 



I 2$2 J 



T{x-Xo±r/KoY 

4$4 ^ y2 



(19) 



A nontrivial feature of expression , which represents 
the entanglement proper, is the combination of two har- 
monics, exp(±zx), multiplied by the Gaussian factors, 
whose maxima are located along two spirals in the plane 
of (r,x): rmax = TKq (x - Xo)- These maxima also de- 
termine the correlation between variables r and x, which 
are an inherent part of the entanglement. The width 
of the maxima gradually spreads out with the growth 
of time (T), proportionally to ^/i^^+T^, as does any 
coherent state evolving in the free space. 

It is relevant here to estimate a dephasing effect of a 
possible mismatch between amplitudes of the colliding 
solitons, A77. It may be estimated through the respective 
change of the relative phase, Ax ^ ijArjAtcoii ~ Arj/c, 
generated by the mismatch during the collision time, 
A^coii ~ l/(?7c) [see Eqs. dJ)-®]. The dephasing is 
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negligible if its size is small in comparison with the per- 
turbation amplitude in expression (|18p . i.e., 



A77/C < e/Ko 



(20) 



Currently available sophisticated experimental methods 
for the creation of solitons in BEC [181 , as well as various 
theoretically elaborated schemes of matter-wave soliton 
lasers [s^l, make it possible to generate nearly identical 
solitons, with a sufficiently small difference between their 
amplitudes. 



C. The initial phase-uniform state 

Instead of initial wave packet (|lip . we can take one uni- 
formly spread over phases (pi and (j)2, which corresponds 
to an experimental situation in which the initial phases 
of solitons are not controlled. Then, modifying expres- 
sion Eq. (fT2|) accordingly, decomposition (fT7|) is replaced 
by its simplified counterpart: 



r+00 

(S/V^) / dKe 



0*°*^ (6,6) = exp 

iK{r-rQ) 



exp 



P 

^2 



{K ~ KoY 



Recombining this with the result of the Born's approxi- 
mation, as per Eq. (I16p . gives rise to the following expres- 
sion for the collision-induced perturbation of the wave 
function: 



S'i (6,6,X) = exp ( ] cosx 



2ie 



i[Koir- ro) - -K^T 



(21) 



Although this result is much simpler than the one given 
by Eq. (fT5|) . in the case when the solitons' phases were 
initially allocated certain values, it is nontrivial too, 
demonstrating the dependence on relative phase x of the 
perturbation generated, in the course of the collision, by 
the phase-independent initial wave function. Actually, 
this feature may be regarded as the simplest manifesta- 
tion of the collision- induced entanglement. The correla- 
tion properties of the entanglement in this approxima- 
tion are obvious, being, as a matter of fact, determined 
by factor cos^. 



IV. CONCLUSION 

In this work, our aim was to present a proof of the 
principle that two initially uncorrelated quantum (actu- 
ally, semi-classical) solitons may get entangled through 
the collision between them. To this end, we have re- 
stricted the analysis to macroscopic coordinates of the 
solitons, viz., the position and phase, which may be justi- 
fied for matter-wave solitons in BEC. In this sense, they 



are regarded as Einstein-Podolsky-Rosen particles with 
the additional internal degree of freedom (the phase) . We 
have presented simple analytic expressions characterizing 
the resulting entangled states, represented by collision- 
induced corrections to initial factorized wave functions 
of the soliton pair. The results were obtained by means 
of the Born's approximation, which is valid for the colli- 
sion between fast soliton. The entanglement and quan- 
tum correlations are described in terms of the relative 
position and phase of the solitons, and the (conserved) 
total momentum and "angular momentum" (the latter is 
the canonical momentum conjugate to the total phase). 
A simple but noteworthy effect is that the collision of 
solitons with the completely uncertain (delocalized) rel- 
ative phase leads to a partial phase localization. The 
predicted results are expected to be directly relevant to 
matter-wave solitons in BEC formed by atoms with at- 
tractive interactions and loaded into a nearly ID trapping 
potential. A realization in terms of optical solitons may 
be possible too, in principle. 

The entanglement predicted by the above analysis is 
weak, as the consideration was limited to the Born ap- 
proximation, i.e., the first approximation of the pertur- 
bation theory. However, under more general conditions, 
the entanglement may reach significant levels (cf. Ref. 
[l3| where the entanglement was analyzed by means of 
an exact quantum solution for NLS solitons interacting at 
zero velocity) . What is then a possible advantage of using 
matter- wave solitons over individual particles, likes pho- 
tons or atoms? The obvious answer is that the solitons 
constitute macroscopic, or at least mesoscopic objects, 
and thus exhibit completely different properties in mea- 
surements. To some extent, they are similar to atomic 
ensembles Q, but, in contrast to the latter, the soli- 
tons represent coherent moving objects. In principle, the 
Bell's inequalities can be tested with entangled solitons, 
and even though such test would not be loophole-free, 
they would evidently exhibit different scaling and level 
of errors, in comparison with previously studied systems. 

As said above, one may attain stronger correlations by 
going beyond the Born-approximation limit. However, 
in that regime the macroscopic (mean-field) description 
may cease to be quantitatively correct - in particular, 
because the deformation of the slowly colliding solitons 
makes the it irrelevant to use the description solely in 
terms of their collective degrees of freedom. Then, it 
may be appropriate to consider quantum solitons using 
the time-dependent GPE [or similar equation(s)] to de- 
scribe the coherent soliton states proper, and the BdG 
(Bogoliubov-de Gennes) equations for quantum fluctua- 
tions around them [sl] , cf. the analysis of the stability of 
matter-wave solitons against quantum fluctuations per- 
formed in Refs. [sO] and [Slj. In particular, care should 
be taken of the coUisional excitation of an IM (intrin- 
sic mode), if the solitons support it. However, such a 
consideration is beyond the scope of this work. 

Similar comments concern another important aspect 
of the studies of the entanglement, viz., the decoherence. 
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In the present part of the decoherence is due to 

the excitation of many-body degrees of freedom (such as 
BdG modes), rather than just quantum fluctuations of 
the sohtons' positions and phases. Nevertheless, in the 
case of the collision between fast solitons we may restrict 
the attention to the collective degrees of freedom of the 
solitons, and look at the respective covariance matrix, 
similarly to the case of Gaussian states. Theii, a pplying 
the entanglement criteria for Gaussian states [53, 0, ISSj 
makes it possible to estimate the level of the entangle- 
ment of the present (ideal) state, as well as the maximum 
allowance of the contribution to the decoherence from the 
neglected many-body modes. 

The estimate of the robustness of the entanglement is 
particularly simple for the solution corresponding to Eq. 
(I^Tj) . because the degrees of freedom corresponding to the 
distance and phase shift between the solitons separate, 
hence the wave function of the two-soliton state can be 
written, in terms of dynamical variables ([7]) and ([5]), as 
[cf. Eqs. (UHl), dm)] 

|M/(p,r,x)) = |*(p))|*(r)) 
x[|0)+C(|l> + |-l»], (22) 

where states of the relative phase (x), |0), | ± 1), cor- 
respond to harmonics with the respective numbers, and 
C EE ej {2\/2K^. The phase part of state ([22]) is ob- 
viously entangled, since its partial transpose has three 
negative eigenvalues, namely, — — Ci^C^ [l^- There- 
fore, the state remains entangled after adding arbitrary 
decoherence (noise) whose density matrix has an opera- 
tor norm smaller than C. This actually means that one 



can add "white noise" with norm C < C to state (|22l) . 
and retain its non-positive partial transpose property. 
As concerns the dephasing due to a possible difference 
in the amplitudes of the colliding solitons, the estimate 
performed above demonstrates that it is negligible under 
condition (PU)) . 

A straightforward extension of the present analysis, 
still within the framework of the description based on the 
mean-field collective degrees of freedom and the Bern's 
approximation, may be elaborated for the description 
of the generation of entanglement by collisions between 
two-dimensional solitons supported by a quasi- ID po- 
tential in the self-attractive BEG, as proposed (in the 
mean- field approximation) in Refs. [sO]- Another exten- 
sion, also possible in the framework of the quasi-particle 
approach, is to consider entanglement induced by colli- 
sions between gap solitons in a self-repulsive BEG loaded 
into an optical-lattice potential. It is known that such 
solitons may be mobile in ID and 2D geometries, with 
a negative effective mass [5l|. The interaction poten- 
tial in the corresponding Hamiltonian is expected to be 
different from the above expression ([S]), namely, being 
^ exp (— |r|) cos(qr) cosx, with some wavenumber g, due 
to the oscillatory shape of tails of the gap solitons. 
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